A stochastic approach to nonlinear unconfined flow subject to multiple random fields

In this study, the KLME approach, a momentequation approach based on the Karhunen–Loeve decomposition developed by Zhang and Lu (Comput Phys 194(2):773–794, 2004), is applied to unconfined flow with multiple random inputs. The log-transformed hydraulic conductivity F, the recharge R, the Dirichlet boundary condition H, and the Neumann boundary condition Q are assumed to be Gaussian random fields with known means and covariance functions. The F, R, H and Q are first decomposed into finite series in terms of Gaussian standard random variables by the Karhunen–Loeve expansion. The hydraulic head h is then represented by a perturbation expansion, and each term in the perturbation expansion is written as the products of unknown coefficients and Gaussian standard random variables obtained from the Karhunen–Loeve expansions. A series of deterministic partial differential equations are derived from the stochastic partial differential equations. The resulting equations for uncorrelated and perfectly correlated cases are developed. The equations can be solved sequentially from low to high order by the finite element method. We examine the accuracy of the KLME approach for the groundwater flow subject to uncorrelated or perfectly correlated random inputs and study the capability of the KLME method for predicting the head variance in the presence of various spatially variable parameters. It is shown that the proposed numerical model gives accurate results at a much smaller computational cost than the Monte Carlo simulation.